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Abstract 

The simulations of field-evaporation processes for surface atoms on W(011) and 
Mo(Oll) surfaces are implemented using first-principles calculations based on the 
real-space finite-difference method. The threshold values of the external electric field 
for evaporation of the surface atoms, which are ~ 6 V/A for tungsten and ~ 5 
V/A for molybdenum, are in agreement with the experimental results. Whereas 
field evaporation has been believed to occur as a result of significant local-field 
enhancement around the evaporating atoms, in this study, the enhancement is not 
observed around the atoms but above them and the strength of the local field is 
much smaller than that expected on the basis of the classical model. 
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1 Introduction 

In the last decade, atomic-scale manipulation on material surfaces using scan- 
ning tunneling microscopy has been the subject of intensive research due to 
its potential for creating artificial surface nanostructures. There have been 
several experiments in which an atom is transferred between a probe tip and 
a surface under the influence of applied voltage pulses (1). This phenomena is 
called field evaporation jj) and is considered to be a thermally activated pro- 
cess in which rate constants can be parameterized according to the Arrhenius 
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formula (3). Since the activation energy for field evaporation varies depending 
on the electric field applied between the probe tip and surface, the contribu- 
tion of the electric field is important during this process. In the past, Tsong 
(0) proposed phenomenological models for evaluating the activation energy for 
the field evaporation of surface atoms. However, for a profound understand- 
ing of this phenomenon from a microscopic point of view, it is mandatory to 
self-consistently calculate both the electronic charge distribution and the in- 
duced electrostatic field. So far, first-principles calculations has been used to 
interpret the filed evaporation of the surface atoms on semiconductor surfaces 
(0;E|). When it comes to the metal surface, although several first-principles 
studies have been implemented (0; Q; H) , all of them employed structureless 
jellium surfaces. 

In this study, first-principles calculations were carried out to explore the ac- 
tivation energies and threshold values of external electric fields for the evap- 
oration of surface atoms from atomically flat W(011) and Mo (Oil) surfaces. 
We determined activation energies for the evaporation of atoms on atomically 
flat (Oil) surfaces under strong electric fields. Our calculated threshold values 
for field evaporation, ~ 6 V/A for tungsten and ~ 5 V/A for molybdenum, 
are in good agreement with experimental data. Whereas field evaporation has 
been believed to occur as a result of significant local-field enhancement (LFE) 
around single metal atoms on atomically flat crystal planes, we observed the 
LFE not around the evaporating atoms but ~ 3 a.u. above them. 

The rest of the paper is organized as follows: In Sec. 2, we briefly describe 
the computational details of our calculations. The results are presented and 
discussed in Sec. 3. We conclude our findings in Sec. 4. 



2 Computational details 



2.1 Method 



Our first-principles simulations are based on the real-space finite-difference 
method Jp; fToh with the incorporation of the timesaving double-grid tech- 
nique fjllh . Compared with the plane- wave approach, the real-space finite- 
difference method is much simpler to implement while maintaining a high 
degree of accuracy. Moreover, the real-space calculations eliminate the serious 
drawbacks of the conventional plane-wave approach, such as its inability to 
describe strictly nonperiodic systems: In the case of simulations under exter- 
nal electric fields using the conventional plane-wave approach, the periodic 
boundary condition gives rise to a saw-tooth potential, which sometimes leads 
to numerical instability during the self-consistent iteration, while the real- 
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space finite- difference method is free from involving the saw-tooth potential 
and can exactly determine the potential induced by the external electric field 
as a boundary condition. The grid spacing is set to be ~ 0.33 a.u. which cor- 
responds to a cutoff energy of <~ 89 Ry in the plane-wave approach and a 
denser grid spacing of ~ 0.11 a.u. is taken in the vicinity of nuclei with the 
augmentation of double-grid points (11). The nine-point finite-difference for- 
mula is used for the derivative arising from the kinetic-energy operator of the 
Kohn-Sham equation (12). The norm-conserving pseudopotentials in- 
corporating a Kleinman-Bylander nonlocal form f|15h are employed to describe 
the electron-ion interaction. Exchange-correlation effects are treated with the 
local-spin-density approximation (|lfih of the density functional theory f|17h . 



2.2 Models 



Figure 1 shows the top view of the (Oil) surfaces employed here. We adopt 
a technique that involves the use of a supercell whose size is set to be L x = 
2^/2^0, L y = 2a and L z = 10a , where L X) L y and L z are the lengths of the 
supercell in the x, y and z directions, respectively. Here, the direction perpen- 
dicular to the surface was chosen as the z direction and the bulk constants a 
are 5.97 a.u. for tungsten and 5.43 a.u. for molybdenum. To completely elim- 
inate unfavorable effects from neighboring cells which are artificially repeated 
in the case of the periodic boundary condition, the nonperiodic boundary con- 
dition of vanishing wave functions out of the supercell is imposed in the z 
direction, while the periodic boundary conditions are employed in the x and 
y directions. The surface Brillouin zone is sampled by the uniform mesh of k 
points and the Fermi level is broadened by the Fermi distribution function. 
The supercell contains three (Oil) atomic layers (i.e., the thin film model) 
and the atoms in the topmost atomic layer are fully relaxed during structural 
optimizations. 



We first determine the minimum-energy atomic configuration of a surface 
atom. A tungsten (molybdenum) surface atom is initially placed at each of 
the sites of A, B, and C in Fig. 1 at an appropriate distance from the W(011) 
[Mo(Oll)] surface, and then the forces on the atoms are relieved. The binding 
energies are listed in Table 1. The most stable site is found to be B on both 
the W(011) and the Mo(Oll) surfaces and hereafter we use this configuration 
as the computational model. 
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3 Results and discussion 



3.1 Activation energy 



We now evaluate the activation energies for the field evaporation of surface 
atoms during the lifting of the surface atoms. One of the potentially important 
applications of the field-evaporation process is the direct determination of the 
binding strengths of surface atoms from the external electric field required for 
their removal. Figure 2 shows activation energies as a function of distance from 
the surface for various external electric fields. These energies become lower 
as the external electric fields are increased. On the basis of these results, the 
threshold values of the external electric fields for evaporation are ~ 6V/A for 
tungsten and ^ 5 V/A for molybdenum, which agree with the experimental 
values of 5.7 V/A for tungsten and 4.6 V/A for molybdenum obtained using 
a field-ion microscope at low temperature. In addition, the threshold values 
of the external electric field are considered to vary depending on the bond 
strengths of the surface atoms, and the rate constants of field evaporation 
depend on the activation energies according to the Arrhenius formula (0); 
the molybdenum surface atom is easily removed compared with the tungsten 
surface atom. This is consistent with the binding energy mentioned in the 
preceding section. 



3.2 Local- field enhancement 



We depict in Fig. 3 the electron distributions p(r, F) — p(r, F = 0) induced by 
the external electric field (18). The overall charge around the surface atoms 
decreases, and the atoms are expected to be positive ions when they evaporate 
(jlflh . We show in Fig. 4 the counter plots of the external electrostatic field 
as the differences between the total electrostatic field in the presence of an 
external electric field and that in the absence of an external electric field 
(jlih . d(V(r, F) — V(r,F = 0))/dz, where V is the sum of the external and 
Hartree potential functions. One can clearly recognize the expulsion of the 
external electronic field from the inside of the surface which is well known in 
the case of macroscopic metallic systems. However the LFE occurs not around 
the evaporating atoms but ~ 3 a.u. above them, and there are no further 
significant differences in external electrostatic field around the evaporating 
atoms. In addition, the strength of the local field is much smaller than that 
expected on the basis of the classical model (16.5 V/A for tungsten and 
13.5 V/A for molybdenum) (Eoh . This situation is consistent with the result 
obtained using a jellium surface (0). 
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4 Conclusion 



We have studied the field-evaporation processes of surface atoms from atomi- 
cally flat W(011) and Mo(OOl) surfaces. The threshold values of the external 
electric field for evaporation of the surface atoms, which are ~ 6 V/A for 
tungsten and ^ 5 V/A for molybdenum, are in agreement with the experi- 
mental results. They depend on the binding energy between the evaporating 
atoms and the surface. Moreover, as the external electric field increases, the 
activation energy for field evaporation becomes lower. While field evaporation 
has been believed to occur as a result of remarkable LFE around the evapo- 
rating atoms, we observe the LFE ^ 3 a.u. above the atoms and the strength 
of the local field is not particularly high compared with those expected on the 
basis of the classical model. 
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Table 1 



Binding energies as the atoms are displaced at sites A-C. 



w 


Mo 


Site 


Binding energy (eV) 


Site 


Binding energy (eV) 


A 


5.09 


A 


3.93 


B 


7.88 


B 


5.65 


C 


6.80 


C 
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Fig. 1. Top view of the two topmost layers of atomically flat (Oil) surface. Large 
and small spheres represent atoms on the top and second layers, respectively. 
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Fig. 2. Potential energies for evaporating atoms on (A) W(011) and (B) Mo(011) 
surfaces as a function of position z of the evaporating atoms. The zero of the position 
z and the energy are chosen to be those of the most stable atomic geometry in an 
absence of the external electric field. 
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Fig. 3. Electronic density shifts p(r, F) — p(r,F = 0) on the (Oil) cross sections 
containing the evaporating atoms on (A) tungsten and (B) molybdenum. Exter- 
nal electric fields F are 5.5 V/A for tungsten and 4.5 V/A for molybdenum. 
The contour spacing is 0.005 electron/supercell and thick curves represent 0.000 
electron/supercell. Solid (dotted) curves represent positive (negative) values. Large 
and small balls indicate the atomic positions on and above the cross section, respec- 
tively. 
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Fig. 4. Differences between the total electrostatic field in the presence of an external 
electric field and that in the absence of an electric field on the (010) cross sections 
containing the evaporating atoms on (A) tungsten and (B) molybdenum. External 
electric fields F are 5.5 V/A for tungsten and 4.5 V/A for molybdenum. The 
contour spacing is 1.0 V/A and and thick curves represent 0.0 V/A. Solid (dotted) 
curves represent positive (negative) values. Symbols have the same meanings as 
those in Fig. 3. 
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